A model reduction approach to numerical inversion for a parabolic partial differential equation
نویسندگان
چکیده
We propose a novel numerical inversion algorithm for the coefficients of parabolic partial differential equations, based on model reduction. The study is motivated by the application of controlled source electromagnetic exploration, where the unknown is the subsurface electrical resistivity and the data are time resolved surface measurements of the magnetic field. The algorithm presented in this paper considers inversion in one and two dimensions. The reduced model is obtained with rational interpolation in the frequency (Laplace) domain and a rational Krylov subspace projection method. It amounts to a nonlinear mapping from the function space of the unknown resistivity to the small dimensional space of the parameters of the reduced model. We use this mapping as a nonlinear preconditioner for the Gauss–Newton iterative solution of the inverse problem. The advantage of the inversion algorithm is twofold. First, the nonlinear preconditioner resolves most of the nonlinearity of the problem. Thus the iterations are less likely to get stuck in local minima and the convergence is fast. Second, the inversion is computationally efficient because it avoids repeated accurate simulations of the time-domain response. We study the stability of the inversion algorithm for various rational Krylov subspaces, and assess its performance with numerical experiments.
منابع مشابه
An inverse problem of identifying the coefficient of semilinear parabolic equation
In this paper, a variational iteration method (VIM), which is a well-known method for solving nonlinear equations, has been employed to solve an inverse parabolic partial differential equation. Inverse problems in partial differential equations can be used to model many real problems in engineering and other physical sciences. The VIM is to construct correction functional using general Lagr...
متن کاملA Nonlinear Parabolic Partial Differential Equation Model for Image Enhancement
Abstract—We present a robust nonlinear parabolic partial differential equation (PDE)-based denoising scheme in this article. Our approach is based on a second-order anisotropic diffusion model that is described first. Then, a consistent and explicit numerical approximation algorithm is constructed for this continuous model by using the finite-difference method. Finally, our restoration experime...
متن کاملParameter determination in a parabolic inverse problem in general dimensions
It is well known that the parabolic partial differential equations in two or more space dimensions with overspecified boundary data, feature in the mathematical modeling of many phenomena. In this article, an inverse problem of determining an unknown time-dependent source term of a parabolic equation in general dimensions is considered. Employing some transformations, we change the inverse prob...
متن کاملNumerical Methods for Fuzzy Linear Partial Differential Equations under new Definition for Derivative
In this paper difference methods to solve "fuzzy partial differential equations" (FPDE) such as fuzzy hyperbolic and fuzzy parabolic equations are considered. The existence of the solution and stability of the method are examined in detail. Finally examples are presented to show that the Hausdorff distance between the exact solution and approximate solution tends to zero.
متن کاملA new approach to using the cubic B-spline functions to solve the Black-Scholes equation
Nowadays, options are common financial derivatives. For this reason, by increase of applications for these financial derivatives, the problem of options pricing is one of the most important economic issues. With the development of stochastic models, the need for randomly computational methods caused the generation of a new field called financial engineering. In the financial engineering the pre...
متن کامل